A Simple Sample Size Formula for Estimating Means 
of Poisson Random Variables * 
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Abstract 

In this paper, we derive an explicit sample size formula based a mixed criterion of absolute 
and relative errors for estimating means of Poisson random variables. 

1 Sample Size Formula 

It is a frequent problem to estimate the mean value of a Poisson random variable based on 
sampling. Specifically, let X be a Poisson random variable with mean K[X] = A > 0, one wishes 
to estimate A as ^ 

n 

where Xi, ■ ■ ■ are i.i.d. random samples of X. Since A is of random nature, it is important 
to control the statistical error of the estimate. For this purpose, we have 



Theorem 1 Let Sa > 0, e,- G (0, 1) and 5 G (0, 1). Then 



Pr 



A- A 



< £a or 



A- A 



provided that 



n> — X 



In I 



Ea {I + er)\n(l + Er) - Er 



(1) 



It should be noted that conventional methods for determining sample sizes are based on normal 
approximation, see [3] and the references therein. In contrast. Theorem 1 offers a rigorous method 
for determining sample sizes. To reduce conservatism, a numerical approach has been developed 
by Chen [1] which permits exact computation of the minimum sample size. 

*The author is currently with Department of Electrical Engineering, Louisiana State University at Baton Rouge, 
LA 70803, USA, and Department of Electrical Engineering, Southern University and A&M College, Baton Rouge, 
LA 70813, USA; Email: chenxinjia@gmail.com 
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2 Proof of Theorem 1 



We need some preliminary results. 

Lemma 1 Let K be a Poisson random variable with mean 9 > 0. Then, Pr{K > r} < e~^ (^)' 
for any real number r > 9 and Pr{/^ < r} < e~^ {t'Y ^'^^ positive real number r < 9. 

Proof. For any real number r > 9, using the Chernoff bound [2], we have 

Fi{K >r} < inf E [e^^^"'^)] = inf V e^^*"^) - 
t>o L J t>o^ 



r 

1=0 



inf e^^*e-^e-^*y ^e-^^' = infe-V 
t>o i\ t>o 

i=0 



where the infimum is achieved at t = In (^) > 0. For this value of t, we have e ^e^^' *^ = e ^ (^)^- 
It follows that Vt{K >r]< e'^ [^)'' for any real number r > 9. 

Similarly, for any real number r < 9, we have Pr{i^ < t} < {~Y ■ 

□ 

In the sequel, we shall introduce the following function 

A 



g{£, A) = e + (A + e) In 



A + e 



Lemma 2 Let A > e > 0. Then, Pr | A < A — e| < exp {n g{—£, A)) and g{—£, A) is monotoni- 
cally increasing with respect to \ ^ (e, oo). 

Proof. Letting K = ^i, 9 = nX and r = n(A — e) and applying Lemma [H for A > e > 0, 
we have 

Pr I A < A - e| = Pr{K < r} < e"^ f y j = exp (n g{-e, A)) , 
where g{—£, A) is monotonically increasing with respect to A € (e, oo) because 



dX \ XJ X 

for A > e > 0. 



□ 



Lemma 3 Let e > 0. Then, Pr|A> A + e| < exp{ng{e,X)) and g{£,X) is monotonically 
increasing with respect to X ^ (0, oo). 
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Proof. Letting K = Xl^Li = nX and r = n(X + e) and applying Lemma [H for A > 0, we 
have ^ 

Pr |a > A + e| = Ft{K > r} < ^ < exp (n g{e, A)) , 

where g{e, A) is monotonically increasing with respect to A G (0, oo) because 

□ 

Lemma 4 g{e, A) > g{—e, A) /or A > e > 0. 

Proof. Since g{£, A) — (7(— e, A) = for e = and 

d[g{e,X)-g{-e,X)] X' , ^ 

5e A2 - £2 ^ 

for A > e > 0, we have 

ff(e,A) -5(-e,A) > 

for any e G (0, A). Since such arguments hold for arbitrary A > 0, we can conclude that 

A) > g{-e,X) 

for A > e > 0. □ 



Lemma 5 Let < e < 1. Then, Pr |a < A(l — e)| < exp (n g{—eX, A)) and g{—eX, A) 



ts mono- 



tonically decreasing with respect to X > 0. 



Proof. Letting K = Y17=i-^i^ 9 = nX and r = nA(l — e) and making use of Lemma [H for 
< e < 1, we have 

Pr |a < A(1 - e)} = Fr{K < r} < < exp (n c/(-eA, A)) , 

where 

g{-eX,X) = [-e-(l-e)ln(l-e)]A, 

which is monotonically decreasing with respect to A > 0, since —e — (1 — e)ln(l — e) < for 
< e < L 

□ 

Lemma 6 Let e > 0. Then, Pr |a > A(l + e)| < exp (n ^(eA, A)) and g{eX, A) is monotonically 
decreasing with respect to X > 0. 
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Proof. Letting K = Y^^=i ^ = ^'^d — nX{l + £) and making use of Lemma[Tl for e > 0, 
we have 



where 



Pr|A>A(l + e)| < exp (n ^(eA, A)) 



5(eA,A) = [e- (l + e)ln(l + e)]A, 



which is monotonically decreasing with respect to A > 0, since e — (1 + e) ln(l + e) < for e > 0. 

□ 



We are now in a position to prove the theorem. It suffices to show 



Pr 



A- A 



>ea& 



A- A 



for n satisfying ([T|). It can shown that ([T]) is equivalent to 

exp(n g{ea,Sa)) < ^. 

We shall consider four cases as follows. 
Case (i): < A < 
Case (ii): A = 
Case (iii): < A < 
Case (iv): A > 



In Case (i), we have Pr{A < A — Ea} = and 



Pr 



A- A 



A- A 



Pr 



A- A 



>ea] 

Pr{A <X-ea} + Pr{A > A + e J 
Pr{A> A + ea}. 



By Lemma 



(2) 



Pr{A > A + ea} < exp(n g{ea, A)) < exp(n g{ea,ea)) < -■ 



Hence, 



Pr I 



A- A 



A- A 



> sA} < 2 



< <5. 



In Case (ii), we have Pr{A < A — = Pr{A = 0} and 



Pr 



A- A 



A- A 



Pr 



A- A 



>ea} 

Pr{A <X-ea} + Pr{A > A + e J 
Pr{A = 0} + Pr{A > A + e4. 
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Noting that ln2 < 1, we can show that —Sa < g{£a^£a) and hence 

Pr{A = 0} = Pr{Xi = 0, i = I,--- ,n} 
= [Pr{X = 0}]'^ 



< exp(n fi((ea,ea)) 



< exp n 5^ 



where the second inequaUty fohows from Lemma Hence, 



Pr<^ A- A 



A- A 



< 



> ErX) < - < 5. 



In Case (in), by Lemma ([2]), Lemma ([3]) and Lemma ([1]), we have 



Pr 



A- A 



>ea 



A- A 



>erA} = Pr{A < A-e4 + Pr{A > A + eJ 
< exp(n5f(-ea,A)) + exp(n A)) 



< exp ( n 5 ( -Ea, — I I + exp I n 5f I — 



< 2 exp g ^e^, ^ < (5. 
In Case (iv), by Lemma ([5]), Lemma ([6]) and Lemma dH), we have 



Pr 



A- A 



A- A 



>e,A} = Prj 



A- A 



> er-A} 

= Pr{A< (l-e^)A} + Pr{A> (l + e^)A} 
< exp(n 5f(— A, A)) + exp(n (7(6^ A, A)) 



Therefore, we have shown Pr 
the proof of Theorem 1. 



A- A 



< exp 5 l^-ea, 

< 2exp^n5rfea, 

> Ea & A - A 



+ exp \ ng [Ea 



< 6. 



< 6 for all cases. This completes 
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